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A COLD FRONT IN A3667: HYDRODYNAMICS AND MAGNETIC FIELD IN 

THE INTRACLUSTER MEDIUM 

A. Vikhlinin*'^, M. Markevitch^, 



(N 

o 
o 

Oh" 

<D 

in 



in 

o\ 
o : 

(N . 
O 

■ 

6 



This conference presentation discusses a Chandra observation of tiie cold front in Abell 3667. We first review our earlier 
results which include a measurement of the front velocity, M ~ 1, using the ratio of exterior and interior gas pressures; 
observations of the hydrodynamic effects expected for a transonic front motion (weak bow shock and gas compression 
near the leading edge of the front); direct observation of the suppressed diffusion across the front, and estimate of the 
magnetic field strength near the front from suppression of the Kelvin-Helmholtz instabilities. 

The new results include using the 2-dimensional brightness distribution inside the cold front (a) to show that the front 
is stable and (b) to map the mass distribution in the gas cloud. This analysis confirms the existence of a dark matter 
subcluster traveling with the front. We also fix an algebraic error in our published calculations for the growth rate of the 
KH instability and discuss an additional effect which could stabilize the front against the small-scale perturbations. These 
updates only strengthen our conclusions regarding the importance of the magnetic fields for the front dynamics. 



A3667 shows clear signs of the on-going merger of two 
big subclusters in its optical, X-ray, and radio images (So- 
dre et al. 1992, Knopp et al. 1996, Rottgering et al. 1997). 
The optical image contains two distinct galaxy concentra- 
tions around the giant elliptical galaxies (we call them A 
and B); the weak lensing observations reveal the associated 
mass concentrations (Joffre et al. 2000). The X-ray emis- 
sion is elongated and co-aligned with the galaxy distribu- 
tion. The most striking feature in the X-ray image is a strong 
surface brightness edge to the South-East of galaxy A, lo- 
cated perpendicularly to the A-B axis. The radio map re- 
veals giant diffuse source of a 0.5-1 Mpc size, located in 
the cluster outskirts and oriented perpendicularly to the A- 
B axis. All available data point to identification of galax- 
ies A and B as the centers of the two merging sublclusters. 
Although it is expected that sublusters should move super- 
sonically (v = 2000-3000 km O"') during the merger, the 
observed line-of-sight velocity difference of galaxies A and 
B is only 120 km s"' (Katgert et al. 1998). This means that 
the merger axis is perpendicular to the line of sight, which 
significantly simplifies interpretation of the data. 

In this conference paper, we first review some earlier re- 
sults from the Chandra observation of A3667 presented in 
Vikhlinin et al. (2001ab). We also present new results on 
mapping of the dark mass of the subcluster traveling with 
the cold front. We also fix an algebraic error in the calcula- 
tions for the growth rate of the KH instability and discuss an 
additional effect which could effectively stabilize the front 
against the small-scale perturbations. 



CHANDRA OBSERVATION OF A3667 

The central part of A3667 was observed by Chandra in 
late 1999. The field of view was centered on the surface 
brightness edge 8' to the South-East of galaxy A. Chan- 
dra angular resolution of 1" corresponds to a proper size 



of 1.46 kpc at the cluster redshift, z = 0.055 (we assume 
Hq = 50 km s"' Mpc"' hereafter). The technical aspects of 
the Chandra data analysis are presented in Vikhlinin et al. 
(2001a). 

The Chandra X-ray image in the 0.7^ keV energy band 
in shown in Fig. |l|. Figure || shows its slightly smoothed ver- 
sion with the main structures marked. The most prominent 
feature is a sharp surface brightness edge whose shape is cir- 
cular with a radius of 410 ± 15 kpc. The edge is very sharp 
— it is unresolved even with Chandra's, angular resolution. 

The temperature map is presented on the bottom panel 
of Fig. |l]. All variations of the temperature in region to 
the North- West of the edge, where the surface brightness is 
high, are significant. To the South-East of the edge the sur- 
face brightness is low and the temperature uncertainties are 
quite large. The very hot spots {T > 12 keV) are statistically 
insignificant. However, it is very clear that the gas temper- 
ature exterior to the edge is much higher than that inside. 
Such structures, observed previously in A2142 (Markevitch 
et al. 2001) were termed the "cold fronts". 

To the South of the cold front, there is a second, 
weaker surface brightness discontinuity clearly visible in the 
smoothed image (Fig. We discuss below that this struc- 
ture can be associated with the bow shock in front of the 
cold cloud moving at a slightly supersonic speed. 

A quantitative study of the cold front properties can be 
made using the observed X-ray surface brightness and tem- 
perature profiles across the front (Fig. H). The brightness 
profile shows a sharp discontinuity. The brightness increases 
by a factor of two within 7-10 kpc of the front position and 
gradually increases by another factor of 2 within 40-80 kpc. 
Such a shape is typical for the projected spherical density 
discontinuity. The deprojection of the surface brightness 
profile gives the gas density n,, = (3.2 ±0.48) x lO'-'cm"^ 
interior to the front, where the error bars conservatively in- 
clude both measurement and all systematic uncertainties. 
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Fig. 3. X-ray surface brightness and temperature profiles across in the ±15° sector across the front. The dashed line in the left panel con'esponds to the 
ROSAT profile without any renormalization. The dotted line corresponds to projection of the spherical density discontinuity. 



Exterior to the front, the gas density is = 0.82 x lO"-' cm"^. 
Across the front, the temperature changes sharply from 
8 to 4.5 keV (Fig. H). The measured temperatures on 
both sides are consistent with constant values within some 
distance from the front. The average temperature within 
275 kpc exterior to the front is Tout = 7.7 ± 0.8 keV, and 
within 125 kpc interior to the front it is ?]„ = 4. 1 ± 0.2 keV. 

Therefore, both gas density and temperature undergo a 
jump across the surface brightness edge. The corresponding 
pressure ratio is pin/pout = 2. 1 ± 0.5 (Table |l]), with a rather 
conservative errorbar. If the cold front were at rest, there 
would be a pressure equilibrium, pi^ = Pout- The higher inner 
pressure means that the front is moving and therefore the 
interior gas "feels" both the thermal and ram pressure of the 
exterior gas. 

HYDRODYNAMICS NEAR THE COLD FRONT 

The observed ram pressure leads to a rather accurate mea- 
surement of the front velocity. Mi = 1.0 ±0.2, i.e. the gas 
cloud moving at the sound speed of the hotter gas. The cloud 
velocity in physical units is v = 1430 ± 290 km s"' . 

Since the cold front moves with a velocity close to the 
speed of sound, the compressibility of the gas should be im- 
portant and therefore two additional structures can be ex- 



Table 1. Gas parameters near the cold front in A3667 



Region 


T 


He 


p = Trie 




IQ-^ cm-3 


10-2keVcm-3 




keV 


Exterior 


7.7 ±0.8 


0.82±0.12 


0.63 ±0.11 


Interior .... 


4.1±0.2 


3.2±0.5 


1.32±0.21 



pected, and indeed observed (Vikhlinin et al. 2001a): a) sur- 
face brightness enhancement due to gas compression near 
the leading edge of the body, and b) a bow shock if the front 
velocity is even slightly supersonic. 

Is the front shape stable? 

In the discussion above, we implicitly assumed that the 
front moves as a solid body. Can the assumption of the 
constant front shape be justified? If the hydrodynamic in- 
stabilities!^ allowed to operate freely, they significantly 
disturb (and even destroy) a gas cloud of radius R by the 
time it travels a distance of several /?'s (Jones et al. 1996). 
Since we observe a very regular shape of the leading edge 
of the cold front, this suggests that the instabilities should be 
somehow suppressed. 

Another indication of the constancy of the front shape can 
be obtained by considering a dependence of gas pressure on 
angle from the direction of motion. All flow lines in the 
immediate vicinity of the front pass through a very small re- 
gion near the leading edge. If the front shape is constant, the 
flow is stationary and therefore the gas pressure should fol- 
low the Bernoulli equation — the pressure decreases as the 
flow accelerates away from the stagnation region. The pres- 
sure inside the front should correspond to the exterior pres- 
sure, otherwise the shape will change and flow will not be 
stationary. Therefore, the assumption of the constant front 
shape leads to an easy and verifiable prediction of the in- 
terior pressure distribution as a function of angle from the 
direction of motion. 

To determine the distribution of gas pressure, we need to 
measure both gas density and temperature as a function of 



In this case the most important type of instability is Rayleigh-Taylor 
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Fig. 1. Chandra X-ray image and temperature map of A3667. Note that the 
surface brightness edge (cold front) is very sharp within the sector 2(/3cr = 
60° , beyond which it is gradually smeared out. 



angle, 9. We find that T = 3.8 ±0.40 keV in the -10° <e< 
10° sector relative to the motion direction, 3.8 ±0.55 keV 
in the -30° <e < -10° sector, and 4.7 ± 0.70 keV for 
10° <6 < 30°.^ Within the uncertainties, these values agree 
with the mean temperature inside the front T = 4. 1 ±0.2 keV 
(Table [l|), therefore we assume that the interior gas temper- 
ature is constant as a function of angle. 

Given the constant temperature, the pressure is simply 
proportional to the gas density which is easily derived from 
deprojection of the surface brightness profiles across the 
front in different sectors from the direction of motion. The 
resulting angular dependence of the gas pressure in shown in 
Fig. H The pressure significantly declines — by a factor of 



^The spectra were collected from a 70 kpc strip just inside the front; 
angles are front North to South through East 
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Fig. 4. Gas pressure inside the front as a function of angle from the leading 
edge. The line is not the fit, but prediction for the M = 1 .05 flow. 



1.5 — at 9 = 30°, where the front starts to widen. Remark- 
ably, the numerical gasdynamic solution (using the VH-1 
code) for the flow about the solid body at Mach numbers 
just above 1 predicts just such a change in pressure (solid 
line in Fig. Thus, the angular dependence of the interior 
gas pressure at least implicitly confirms that the front moves 
with a constant shape. 

Note that if we can assume a priori that the front shape 
is stable and the gas temperature inside is uniform, the de- 
pendence p(9) is a convenient diagnostics of the front speed. 
Indeed, p{9) is rather sensitive to the flow speed; for exam- 
ple, the shape of the measured angular dependence of pres- 
sure shown in Fig. ^ can be reproduced only for 0.7 < M < 
1.5. Under the assumption of constant temperature, p(9) 
is equivalent to p{9) / p(Q), which is easily derived from the 
imaging data only. 
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Fig. 5. Confidence intervals(68% and 95%) for location of the centroid of 
the King profile d (distance from the center of curvature of the cold front), 
and total mass within r = 410 kpc. 



Dark halo associated with the cold front 

Since the shape of the cold front remains constant, we can 
naturally assume that the interior gas is in hydrostatic equi- 
librium in the gravitational field in its reference frame. In 
this case, the hydrostatic equilibrium equation, which under 
assumption of constant temperature reduces to 



p= Pq expi-finipip/T), 



(1) 



allows to infer the distribution of mass of the subcluster trav- 
eling with the front. We assume that the subcluster dark 
matter density profile is given by the King profile, p,„ = 
Pmoi l + r^ / r^Y^I'^ . The density distribution of the cold gas is 
then given by equation (|l]) where iy9 = 47rG/3o ( 1 - asinh jc/jc) , 
denoting x = r/rc- Integration of the plasma emissivity, 
£ oc p^, along the line of sight gives the model distribution 
of the surface brightness. 

The maximum likelihood fit of this model to the ob- 
served 2-dimensional surface brightness distribution within 
approximately 200 kpc of the front gives the parameters of 
the King profile — the distance, d, of the dark halo cen- 
troid from the center curvature of the front, its core-radius 
rc, and central density which we parameterize as the total 
mass within a sphere of the 410 kpc radius (the radius of the 
front). The fit is restricted to the region near the leading 
edge of the body due to the following reasons: 1) this re- 
gion is most interesting from the point of view of the front 
dynamics; 2) we can reasonably assume that the gas dis- 
tribution there is axially symmetric; 3) this region is small 
enough for the King profile to have sufficient freedom for 
approximating any reasonable distribution of mass. 

The results of the maximum likelihood fitting are shown 
in Fig. ^ We find that the center of mass is located at 




equipotential surfaces 



Fig. 6. 



c/ = 315 ±7 kpc from the center of curvature of the cold 
front (i. e. only at 95 kpc from the leading edge), and the 
total mass within 410 kpc is (3.2 ± 0.8) x IO'^Mq. The al- 
lowed values of the core-radius are in the range 50 kpc < 

< 200 kpc. Interestingly, the derived values of mass and 
core-radius are very reasonable for the central region of a 3- 
4 keV cluster. The total mass within a sphere of the 100 kpc 
radius (the distance to the front) is 4.6 x lO'^M©, with ex- 
ceeds the gas mass within the same radius by a factor of 15. 

The distance from the subcluster center of mass to the 
leading edge of the front is only 25% of the front radius. 
This means, roughly speaking, that the dark matter halo 
moves in front of the cold gas and pulls it along. Note that 
such a configuration is required for the global stability of the 
front, as illustrated in Fig. ^. Indeed, in the absence of any 
gravitational field, the gas pressure inside the front should be 
constant, while the exterior pressure should decrease from 
the leading edge to the side according to the Bernoulli equa- 
tion. Therefore, the front shape could not be stable. The 
gravitational potential of the dark matter halo allows a sta- 
ble non-constant pressure distribution of the interior gas to 
emerge. It also provides stability against the global pertur- 
bations of the front shape. Indeed, let the front bend so that 
its curvature increases. The exterior pressure at the outskirts 
will decrease because the flow speed increases. At the same 
time, the interior pressure will increase since the front will 
be at a deeper gravitational potential level. The resulting 
pressure difference will restore the original front shape. 

The presence of the massive dark halo is necessary also 
for suppression of the small-scale Rayleigh-Taylor instabil- 
ities. Due to ram pressure of the ambient gas, the cold cloud 
slows down, which leads to the effective gravity force di- 
rected from the inside of the cold front outward. Since the 
ambient gas density is lower, this configuration is unstable. 
Let us make numerical estimates. The drag force on the 
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Fig. 7. (a) Illustration of the formation of the magnetic layer near the front surface. The initially tangled magnetic lines in the ambient hot gas (red) are 
stretched along the surface because of tangential motion of the gas. The magnetic lines inside the front are stretched because, in the absence of complete 
magnetic isolation, the cool gas experiences stripping. This process can form a narrow layer in which the magnetic field is parallel to the front surface. 
Such a layer would stop the transport processes across the front, as well as further stripping of the cool gas. (h) The interface between the cool and hot gas 
is subject to the Kelvin-Helmholtz instability. The magnetic layer can suppress this instability in the region where the tangential velocity is smaller that a 
critical value Vcr. The velocity field shown for illustration corresponds to the flow of incompressible fluid around a sphere. 



cloud is F(j = C X PomV^A /2, where pout is the density of the 
ambient gas, A is the cloud cross-section area, and C k, 0.4 
is the drag coefficient for a cylinder with a rounded head at 
transonic speed. The cloud acceleration is then 



0.2 Pout V^Trr^ nicPoutV^ 



4/37rrVi„ 



(2) 



Substitution of numerical values yields « 
8 X 10~'°cm s"^. At the same time, the gravitational 
acceleration near the front due to the subcluster mass is 



GM{< lOOkpc) 
(100kpc)2 



:6.4x 10"'' cm s-^ 



(3) 



acting in the opposite direction, and therefore the accelera- 
tion due to drag is negligible. Note that in the absence of 
dark matter, the gravitational acceleration (due to the gas 
only) would be a factor of 15 smaller, which is insufficient 
to compensate the drag acceleration. Therefore a purely gas 
could would be Rayleigh-Taylor unstable, as is indeed seen 
in numerical experiments (Jones et al. 1996). 

MAGNETIC STRUCTURE NEAR THE COLD FRONT 

Magnetic fields can profoundly affect the properties of the 
intergalactic medium (IGM) in clusters. The current mea- 
surements of intracluster magnetic fields are based on Fara- 
day rotation in radio sources seen through the IGM (e.g.. 



Kim, Kronberg & Tribble 1991), or on combined radio and 
hard X-ray data on cluster radio halos under the assumption 
that the X-rays are produced by inverse Compton scatter- 
ing of the microwave background (e.g., Fusco-Femiano et 
al. 2000). Both these methods indicate the magnetic field 
strength on the microgauss level, with considerable uncer- 
tainty. The properties of the cold front in A3667 allow to 
determine the magnetic field strength by a completely new 
method, from its effect on dynamics of the intracluster gas. 

The idea is as follows. We have seen that the cold front is 
very sharp and has a smooth shape; therefore, it must be 
mechanically stable. At the same time, it can be shown 
that given the observed velocity of the ambient gas flow, 
the sharp front must be rapidly destroyed by the Kelvin- 
Helmholtz instability. An examination of the possible mech- 
anisms for the suppression of this instability points to the 
magnetic field surface tension as the most probable candi- 
date. 

The required magnetic field configuration can form 
through a scenario schematically shown in Fig. 0. An ini- 
tially random magnetic field is frozen in the intracluster gas. 
Because of the tangential plasma motions around the front, 
the magnetic loops are stretched along the surface. The mag- 
netic reconnection can then produce a layer in which the 
field lines are parallel to the interface between the cold and 
hot gases and to the flow. This configuration of the magnetic 
field also provides surface tension that can suppress the KH 
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instability. 

The flow speed increases as a function of distance from 
the cold front. If the magnetic field is not too strong, sooner 
or later it will be unable to stabilize the KH perturbations, 
and at this point the instabilities will start to grow; the loca- 
tion of this point is marked as (^^cr is Fig. ^ We expect that 
within the sector ztipd- is sharp, while beyond it, the cloud 
boundary is smeared by turbulence. One readily notes the 
analogy with the X-ray image of A3667, where the front is 
sharp within ip « 30° and quickly smears beyond this radius 
(Fig. Assuming that at (/j = 30°, the magnetic field ten- 
sion is just enough to provide mechanical stability, we can 
determine the magnetic field strength. 

Let us consider in detail the line of argument. For sim- 
plicity, we assume that the front velocity is M = 1 exactly. 

Hydrodynamic KH instability of the cold front 

The interface between two fluids in most cases is sub- 
ject to two types of hydrodynamic instability: (i) Rayleigh- 
Taylor instability, which arises if the effective gravity is di- 
rected from the heavier fluid into the lighter one; and (ii) 
Kelvin-Helmholtz instability, which arises from tangential 
motion of the fluids. As was discussed above, the dark 
matter halo moving with the front suppresses the Rayleigh- 
Taylor instability. Therefore, we consider below only the 
KH instability. 

A numerical solution of the M = 1 flow about the cold 
front gives the following angular dependence of the velocity 
of the ambient gas near the front: 

M((/7) ~ 1.1 sin(^. (4) 

At (y5 = 30° the Mach number reaches M = 0.55. Clearly, this 
is a high velocity and therefore one expects a fast growth 
of the KH perturbations leading to formation of a turbulent 
layer Below, we will use the following widely known prop- 
erties of the KH instability: 1) The fastest-growing perturba- 
tions are parallel to the flow; 2) The short-wavelength modes 
grow faster than the long-wavelength modes; 3) The pertur- 
bations grow exponentially, and therefore a nonlinear stage 
is reached rapidly. When the nonlinear stage is reached, the 
amplitude is approximately equal to the wavelength of the 
initially linear perturbation; subsequently, the perturbations 
with a given wavelength merge into large vortices. Thus, if 
one can show that the perturbations with wavelength A are 
unstable, the turbulent layer of width at least A is expected. 

Dispersion equation 

The stability analysis is conveniently performed via con- 
sideration of the small perturbation of the interface in the 
form of the planar waves and solving the resulting disper- 
sion equation. In our case the dispersion equation can be 
written as (Miles 1958): 

-±- -^/-g +^=0 (5) 



Hot gas 

Ph, Th, Ch 




V 

M = V/Ch 


Cold gas 

Pc, Tc, Cc 







Fig. 8. Hot gas flows with a velocity v" relative to the colder, denser gas. The 
interface is subject to the planar wave perturbations, exp(!(fcc — uit)), with 
the wave vector parallel to the flow. For the given wave vector, ijj can be 
found from the dispersion equations: for pure hydrodynamic and ( p^ ) 
in MHD. 

where M is the Mach number of the hot gas, and c^ and 
ci, sound speeds in the cold and hot gases, respectively (see 
Fig.l). 

Application to A3667 

For all relevant Mach numbers, equation (^ has two com- 
plex solutions for oj, one of them corresponding to the grow- 
ing mode. The instability growth time is r = (Imcj)"^ Al- 
though perturbations with any wavelength are formally un- 
stable, only for some of them the growth time is sufficiently 
short compared to the lifetime of the front. The relevant 
time scale with which r should be compared is the cluster 
core passage time, fd-oss = L/M^c^, where L is the cluster 
size. The value of exp(t„oss/T) is the perturbation growth 
factor over the core passage time; if fcross/T ^ 1 (we will 
use a conservative fiducial threshold of 10), the perturbation 
is effectively unstable because it has enough time to become 
non-linear 

For the observed gas parameters and = 1 the solution 
of (^ is approximately given by fcross/T — 12ML/ A, where 
A = 27r/fc is the perturbation wavelength. Using the distri- 
bution of the Mach number near the front (eq. we have^ 

'cross o o ^ • 

~3.3--smo9. (6) 

T A 

This equation, applied for a typical cluster size of L = 1 Mpc, 
implies that perturbations with A = 10 kpc are unstable 
(fcross/T = 10) at (f > 1.7°, that is, essentially everywhere; 
at (p= 10°, the modes with A < 57 kpc are unstable, and at 
ip = 30°, modes with A < 150 kpc. 

There is another effect which limits the lifetime of per- 
turbations. The linear perturbations have non-zero group 
velocity Vdr = dRe lu /dk, and therefore drift along the flow 
as they grow. In principle, in some cases the perturbation 
can move to the side of the front before it grows non-linear. 
A detailed analysis of solutions of eq. (^ shows that uj is 
always proportional to k and therefore the group velocity 

^We note that the analogous equation (4) in Vikhlinin et al. 2001b con- 
tained a factor of (27r)^ smaller numerical factor due to an algebraic eiTor. 
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equals the phase velocity, Vdr = Re(a;)/^. Furthermore, in 
the entire range of the flow velocities near the cold front in 
A3667, the imaginary part of uj is almost proportional to its 
real part: 

^^0.70. (7) 

This equation implies a simple relation between the pertur- 
bation growth factor and its displacement, A(p: 

Growth factor = 



exp / Imujdt ~ exp 1.43 / Reujdt 



:exp( lA3k I ^Y~dt\ =e\pflA3k f dl^ 



where R is the radius of the front. Assuming, as before, that 
reaching a growth factor of e'" is required for a perturbation 
to become non-linear, we obtain the following condition for 
effective instability for /? = 410 kpc: 



A < 190 kpc 



Aip 
30°' 



(9) 



The longer wavelength perturbations drift to the side of the 
front before they grow non-linear Perturbations with the in- 
teresting A ~ 30 kpc should grow strongly non-linear before 
they leave the observed smooth region of the front. 

To summarize, the perturbations with A < 30 kpc are KH- 
unstable over most of the front surface, and therefore a for- 
mation of the turbulent layer of at least this width is ex- 
pected. Since the observed front width is < 5 kpc (Vikhlinin 
et al. 2001a), the KH instability must be suppressed. 

The commonly known mechanisms for suppression of the 
KH instability are a) gravity and b) surface tension at the 
interface. Let us consider them in application to A3667. 

The role of gravity in suppressing the KH instability 

As we showed above, there is a massive subcluster travel- 
ing with the front and its gravity has some stabilizing effect 
on the front. The gravitational force suppresses the KH in- 
stability if it acts in the direction from the lighter fluid into 
the heavier one, and if the following stability condition is 
satisfied: 

2 PcPh 



Pt 



-Pi 



It is clear from this equation that only the long wavelength 
modes are suppressed by gravity, while the shorter wave- 
length perturbations can grow almost freely. Noting that 
Pc ^ Ph, equation (nOh can be rewritten as 



> ,,2^/1 2 2 



Th 



'-pmp/l Pmp/l 



(11) 



The gravitational acceleration in the front rest frame was de- 
rived above. Using its numerical value, we obtain the fol- 
lowing stability condition: 



A 



1 kpc 



> 3500 = 4200 sin^ if 



(12) 



This implies that the perturbations with A = lOkpc are sta- 
bilized by gravity only within the sector cp < 3°; at ip = 5° 
stability is achieved only for A > 32 kpc, and at ip = 30° — 
for A> 1050 kpc. 

Therefore, gravity cannot suppress a formation of a turbu- 
lent layer with a width of 10-20 kpc. This leaves the surface 
tension of the magnetic field as the most likely mechanism 
for suppression of the KH instability. 

Suppression of the KH instability by magnetic field 

The KH instability can be suppressed by the surface ten- 
sion of a magnetic field, if the latter is parallel to the in- 
terface and to the direction of the flow. We have discussed 
above that even if the magnetic field is initially highly tan- 
gled, the tangential gas motion near the front can create the 
required magnetic configuration (see. Fig. Q). 

The analysis of KH instability in MHD is greatly simpli- 
fied if the plasma can be considered incompressible. Fortu- 
nately, this can be done in our case because the gas exterior 
to the front flows with a relatively low local Mach number 
M < 0.5, and the growing modes of the KH instability in the 
interior cool gas generally have low phase speed. The in- 
compressibility assumption can be also justified by the fact 
that the growth time of the hydrodynamic instability given 
by eq. (Q) is very close to that in the incompressible limit, 
i. e. for Mc;, -> v, c^/cj ph/pc, Cc -> oo. 

The dispersion equation for small perturbations in a per- 
fectly conducting, incompressible plasma can be written as 
(Syrovatskij 1953) 



■g2 d2 

p,,{u^-kvf + pcu:^ = k'[^ + ^ 



(13) 



where Z?/, and Be are the magnetic field strengths in the hot 
and cold gas, respectively (we assume that the field lines are 
elongated with the flow), and ph and pc are the gas densities. 
The solutions of this equation are real (the interface is stable) 
if 

(14) 



PhPc 



Ph + Pc 



(10) Application to A3667 



If the magnetic pressure is small compared to the thermal 
pressure, p (as it turns out to be the case), or is the same 
fraction of p on both sides of the interface, the stability con- 
dition ( p4| ) can be rewritten in terms of the thermal pressure 
and temperature of the gas: 



B| -B? 1 



(15) 
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Therefore, for the observed temperatures and flow veloci- 
ties the stability of the cold front within the sector tp < 30°, 
where M ^ 0.55, requires that {BI + BI)/2,tt > 0.11 p. If 
smearing of the front beyond this sector is interpreted as the 
onset of the instability, this becomes an estimate of the total 
magnetic pressure: 

/'mag.A+Pmagx = 0.17pgas (16) 

The formal statistical uncertainty in the derived value of 
the magnetic field strength is mostly due to uncertainty in the 
angle, ip„, where the interface becomes unstable, and hence 
in the local Mach number of the hot gas. For ipa in the range 
30° ±10° (which appears to be a conservative interval, see 
Fig. |l]), we find from eq. (|l4|) and (Q) that the magnetic field 
strengths are in the interval 

0.09p<-^ '-<0.23p. (17) 

Stt 

The maximum of the magnetic field strengths in the cold and 
hot gases, B = max(B/,,Bc), is in the interval (47r;?n,ag)'''^ < 
B < (STrpmag)'/^, where p^ag = (B^ + B^)/87r is constrained 
by (|l7|). Using the value of gas pressure inside the dense 
cloud, we find the corresponding uncertainty interval 7 |jG < 
S< 16nG. 

Thus, the front sharpness near the axis and its gradual 
smearing at large angles are most likely explained by the ex- 
istence of a layer with a ^ 10 \iG magnetic field parallel to 
the front. Equation ( [l6| ) shows that the magnetic pressure in 
the layer is a small fraction of the gas pressure. Everywhere 
else it should be even smaller, because the magnetic field in 
the layer is amplified by stretching of the field lines. This 
conclusion is robust: e.g., if the magnetic pressure were of 
order of the thermal pressure, eq. ( p3| ) implies that the front 
would be stable for Mach numbers up to M ^ 1.9, i.e. over 
the entire surface of the cool gas cloud. 

Faraday rotation measurements of the magnetic field out- 
side the cluster cooling flows provide the values of B from 
~ 1 |jG (Kim et al. 1991) to - 10 |jG (Clarke et al. 2000), 
with considerable uncertainty of the absolute field strength 
due to unknown degree of entanglement. Near the cold 
front, the field is straightened and amplified, so our estimate 
likely is an upper limit on the absolute field strength in the 
cluster core, and thus is consistent with the Faraday rotation 
measurements in the bulk of IGM. 

Note that to suppress the instability, the magnetic field 
does not have to be completely ordered. It is sufficient to 
arrange for uniformity on the scales of order 100 kpc. The 
perturbations with A < 10 kpc will then stabilized, while the 
growth time for perturbations withzA > 100 kpc is longer 
than the front lifetime anyway. 

The stretching of magnetic field lines by tangential 
plasma motions is not the only mechanism for creation of 
a layer with a highly ordered magnetic field. For exam- 
ple, Frank et al. (1996) and Jones et al. (1997) showed 



that even a weak seed magnetic field which is incapable of 
suppressing the KH instability, is significantly amplified in 
the nonlinear vortices, and after active field line reconnec- 
tions, a thin layer is formed. In this layer, the field lines 
are aligned with the flow and the field strength increases 
so that interface becomes stable. Therefore, the formation 
of a magnetic layer considered above seems natural; once 
formed, the layer acts as "magnetic isolation" which stops 
the transport processes across the front, and also as "mag- 
netic lubrication", which provides for nonviscous, laminar 
flow around the cold front. 

In conclusion, we note that the structures similar to the 
cold fronts in A3667 and A2142 have now been detected in 
many other clusters; this is a rather common phenomenon 
(see Markevitch et al. 2002 for a review). 

We thank P. Mazzotta for pointing out an algebraic error 
in derivation of equation (6). 
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